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In 2D XY ferromagnet the dipole force induces a strong interaction between 
spin-waves in the long-wavelength limit. The major effect of this interaction is the 
transformation of a propagating spin-wave into a diffusion mode. We study the 
anomalous dynamics of such diffusion modes. We find that the Janssen-De Dominies 
functional, which governs this dynamics, approaches the non-Gaussian fixed-point. 
A spin-wave propagates by an anomalous anisotropic diffusion with the dispersion 
relation: iuj^k^'" and iuj^k^"^ , where = 47/27 and A^^ = 47/36. The low- 
frequency response to the external magnetic field is found. 
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I. INTRODUCTION 



The notion of the ground state and independent elementary excitations hes in the back- 
ground of the condensed matter theory. The excitations, such as electrons, holes, phonons, 
excitons, spin-waves etc. normally have a propagating, wave-like nature. In a homogeneous 
medium they are characterized by their momentum or quasi-momentum p and their energy 
u). Each kind of excitations has its specific dispersion relation or spectrum oj = e(p). A wide 
scope of physical problems can be solved assuming the excitation to be independent, or con- 
sidering their interaction as a weak perturbationB. However, in recent years more and more 
problems have occurred to go beyond the simple picture of non-interacting or weakly inter- 
acting excitations. In his pioneering work A.B.Migdali has indicated that electron-phonon 
interaction is not weak in a narrow range of energy, leading to a strong renormalization of 
the Fermi- velocity and even to an instability of the Fermi-surface. 

Recently a growing number of physical systems revealed excitations which bare spec- 
trum, obtained from linearized equations of motion, is strongly distorted by interaction 
with vacuum and thermal fluctuations. A few recent examples are: the so-called marginal 
Fermi-liquid, in which electrons interacts via transverse magnetic field fiuctuationsi, and 
z/ = 1/2 state of the Fractional Quantum Hall effect in which initial electrons transform into 
quite different fermionsi. 

On the other hand, many years ago, the mode-mode interaction has been recognized 



be considered as hydrodynamic modes in the long- wavelength limit, as well as propagating 
waves, such as sounds, spin waves etc. Their interaction has been proven to be substantial 
not only in the critical region, but also for the hydrodynamics of liquid crystals!'! and for 
the CDW phason modes interacting with impuritie^l. The excitation spectrum of these 
systems is reconstructed by effects of strong interaction. 

Here we present a new solvable and experimentally feasible situation where the strong 
interaction between spin-waves leads to the replacement of the propagating spin-wave by a 




Particles, heat and spin diffusion must 
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diffusion mode and to the appearance of a new soft-mode in a certain range of momentum. 
This is the two-dimensional XY ferromagnet with dipolar interaction between spins. 

The spin-diffusion mode appears naturally in the paramagnetic phase and in the vicinity 
of the Curie pointS. We consider, on the other hand, a low temperature ordered phase, 
where no diffusion is expected but rather a propagating and weakly dissipating spin-wave 
mode. 

In a 2D XY ferromagnet at low temperatures the dipolar interaction is relevant in the 
long wavelength limit, even despite of the low density of spin- waves. It was shown by one 
of the authorJli that dipolar force induces an anomalous anisotropic scaling of spin-spin 
correlations in the ordered phase. In this article we find an analogous dynamical scaling. 

This paper is organized as follows. In the next section we define the model and describe 
the spin-wave spectrum. In section III we discuss the dynamics of XY-magnet with the dipole 
interaction and formulate the perturbation expansion for the model, using the Janssen-De 
Dominic^i technique. Section IV is devoted to the solution of the Dyson equation. There 
we find the self-induced dissipation of the spin-waves. In section V the renormalization of 
the diffusion mode is considered and the anomalous anisotropic exponents are found. In 
section VI the dynamical susceptibility is found. In conclusion we discuss prospects of the 
experimental observation of the anomalous modes. In the Appendix A we use the Ward- 
Takahashi identities to prove that the vertex corrections are small. Details of the Dyson 
equation solution can be found in the Appendixes B and C. 

A brief report on main results of this article has been published earlieilli. 



II. THE HAMILTONIAN AND THE SPIN- WAVE SPECTRUM. 

Spin-waves are fundamental excitations of the exchange magnet with spontaneously bro- 
ken continuous symmetry. According to the so-called Adler principle^l, the interaction 
between spin waves vanishes in the long- wavelength limit. At low temperatures the equi- 
librium density of spin waves is relatively small, and a long wavelength non-equilibrium 
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spin wave, excited by an external source, decays into other spin- waves or scatters on an 
equilibrium spin wave slowly. Thus, the dynamical properties of the exchange magnet are 
determined by the well-defined spin-wave mode. It also means that the imaginary part of 
the poles of the dynamical response function becomes much smaller than their real part as 
the wavelength grows to infinity. 

The 2D exchange magnet with the easy-plane anisotropy is described by a following 
classical Hamiltonian: 



H 



ExA\ 



Jd'x -(VS)2 + -^2 _^^^^|d\SH, (1) 



where J is the exchange coupling constant, A is the strength of the easy-plane anisotropy 
(A > 0)0, g is the dipole interaction coupling constant, H is the external magnetic field, go 
is gyromagnetic ratio, /^s is the Bohr magneton and the field S(x) represents the local spin 
of the magnet and can be normalized by a constraint 

S^(x) = 1. (2) 

The magnetic dipolar energy Hdip is represented by the sum: 

where £> is a unit vector pointing from Xj to Xj, Sj = S(xi), and g = 27i{gGfiBSa~'^y for 
the square lattice; for other lattices g = 2n{gGftBScr)'^, where a is the inverse area of a 
plaquett of the lattice. The magnetic dipole energy can be separated into a short-range and 
a long-range part in the standard way0. The short-range part renormalizes the single-ion 
spin anisotropy and favors in-plane spin orientation. The long-range part of the magnetic 
dipole energy is conveniently expressed in terms of the Fourier-transform Sk of the local 
magnetization field S(x), 

rr 9 (Skk)(S„kk) , . 

^dip = -2^ l^j ' [V 

Thus, the total magnetic Hamiltonian is: 



H = HexA + Hdip. (5) 

The dipole forces are crucial in the long- wavelength limit. Under the scaling transfor- 
mation the short-range exchange interaction scales like L'^~^, whereas the dipole force scales 
like L'^'^^^, where c? = 2, 3 is the spatial dimension of a magnet and L is the scale. The dipole 
energy is characterized by the dipole constant g and the exchange energy is characterized by 
the exchange constant J. Normally in ferromagnets ga'^~^ is much smaller than J (a is the 
lattice constant). However, as the scaling transformation shows, beyond the characteristic 
scale Ld ~ {J / g) in 2D and Ld ~ \J J / 9 in 3D the dipole interaction dominates the energy 
of a spin- wave. 

In the 2D XY ferromagnet the dipole force stabilizes the long-range orderlll. According 
to the Landau-Peierls-Hohenberg-Mermin- Wagner theorem (see e.g.El ) in the absence of 
a long-range interaction the 2D magnet with broken continuous symmetry exhibits the 
algebraic decay of the spin correlations instead of the long-range order, and the infinite 
susceptibility at zero magnetic field. Thus, the dipole force plays a special role in 2D XY 
ferromagnet. In addition, the dipolar force is crucial for the spin statics and dynamics in 
the ordered phase. In contrast to the 3D ferromagnet, where the spin waves are almost 
free, in 2D the interaction between spin waves induced by dipole force is dominant in the 
long- wavelength limit. 

The Hamiltonian @ has two different scales: the anisotropy scale La = \fjj^ and the 
dipole length L/j = J/ g. We assume that the anisotropy A is large compared to the dipolar 
energy {La < Ld). We direct the y-axis along the net magnetization of the magnet and the 
z-axis perpendicular to the plane. The unit vector field S can be represented by two scalar 
fields 0(x, t) and 7r(x, t) 

S = (—V 1 — TT^ sin 0; Vl — TT^ cos 0; vr^ , (6) 

where both vr and are small due to the fact that the dipole force stabilizes the long-range 
order. 

With the precision to the fourth power of and vr the Hamiltonian @ is: 
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k,ijj 



~k,— tij 



(7) 



Here (j)^^^ and TTk^^j are the Fourier-transforms of the fields 0(x, t) and 7r(x, t) respectively. 
We expanded the in-plane magnetization components cos(f) and sin (p up to the fourth power 
in small spin fluctuations (p. We take the uniform magnetic field H to be directed along the 
y-axis and h = gcfJ^BSH. The Fourier-transformed quantities are defined by 



(x,t)= 



where an abbreviated notation 



d'k du 



(8) 



is used and 



■kJu J J (27r)2 27r 
denotes the Fourier-transformation of ^-4^. 



The interaction between the spin waves is described by the non-quadratic part of the 
Hamiltonian (|^: 



mt 



/(kl, k2, k3)0ki,C^i0k2,<^2 0k; 



3,^3 



ki,k2,k3 J uji,ui2,ui3 

x5(ki + k2 + ks)6{uji +UJ2 + cus) 



+ / / ^(ki,k2,k3,k4)0k 

kl ,k2 ,k3 ,k4 LJi ,u:2 ,oJ3 ,1^4 



x5(ki + k2 + ka + k4)5(u;i +UJ2 + UJ3 + W4) 



(9) 



where the three-leg bare vertex /(ki, k2, ka) is: 



/(ki,k2,k3) = |E7^ 



and the four-leg vertex is defined as: 



24 |k, + k,| 



(10) 



(11) 



The vertices ( ]T(]| ) and (|TT]) decrease as the momentum k goes to zero. Besides of that, 
these vertices are singular at |k| 0. Nevertheless, the interaction between spin waves 



asymptotically vanishes in the long wavelength limit for a 3D mag netil. No renormalization 
of the bare correlator 

T 



ir(k) = (55,(k)<55,.(-k)) 



c) 



h.2 

|k| 



(12) 



appears in this limit (we put h = 0). 

In 2D the situation changes drastically: the interaction grows with the wavelength, 
resulting in strong renormalization of critical exponents. To show the difference we calculate 
the upper marginal dimension of the Hamiltonian (|^). Let us consider an arbitrary diagram 
from the perturbation expansion of some correlator. In order to add an internal line to such 
a diagram, we need to add three bare correlators K given by eq. (0), two vertices / and 
one integration over dk^S^-^hy. Fromeq. (0) we see that for small momenta, ~ ky^^ 
and K ~ A;"^. Equation (jT^) gives f ^ k^ ^ ^^y"^ ■ Hence, if we require that the one-hne 
insertion be dimensionless, we obtain: 3 ■ (—2) + 2 ■ 3/2 + 3/2 + (/} — 1) = or D = 5/2. It 
means that the theory is renormalizable in < 5/2. It turns out that in the static 2D case 
the critical exponents can be found exactly. 

Following the worklll, we rescale the field ^ — > ^ j \fJg and rewrite the Hamiltonian (^ 
in a slightly different form: 



akl(f)k(f)_k + a 



-1 k.. 



+ 2w- 



+ w'^alk^ 



J -k 



J k 



(13) 



J -k/ 



where a = ^J J/g, = \J g/ and we have taken into account that for small k one can 
substitute |k| —^ky (see (|12D). Requiring that the Hamiltonian (|13D does not change under 
scale transformation ky Iky, k^ l^^k^, 0k ^^'''4>k, cl l^°a, w l^°w, we find 
the bare exponents A° = 3/2, A° = 0, A° = —9/4, A° = 1/4 (where we have taken into 
account that ^ has one more integration). Now, according to the standard procedure^, we 
introduce the renormalized field and charges: 0^ = ^</,0jj, a = ZattR, w = Z^fi^^'^wji, where fx 
is the scale at which the dimensionless wr is observed. The second term in the Hamiltonian 
(P!B|) is non-analytic at small momenta. It means that there can be no corrections to this 
term in the regular perturbation expansion. Furthermore, from Appendix A one can see. 



that in this case, due to symmetry, there will be no correction to the third term either. 
Hence, we can establish the relations between the renormalization coefficients: 

^a'Z^ = 1 (14) 

Zfz^fj.'/' = 1. (15) 

It can be proven that the renormalization group flow has a non-zero fixed point. It means 
that wr =const and hence ~ (A/zi)^/"^, where A is the scale at which the coupling 
constant equals w. Introducing the critical exponents ~ A^^**, Za ~ A'^", we find from 
(14) and ([T5|) that = —1/12 and = —1/6. Demanding that the first and the second 



terms in the Hamiltonian (|T^) have the same dimensions, we find A^^ =4/3. Finally, from 
the Callan-Symanzik equationS, the long range limit of the two point correlation function 
can be found: 

Ki^) = {x' + \yn''/'f{x/\yn, 

where f{x) is an arbitrary finite function. The dynamical properties of the magnet, however, 
are still to be found. 

III. THE DYNAMICS OF XY MAGNET 

At T = the classical magnet with the Hamiltonian (^, obeying the constraint (^, 
follows the Landau-Lifshitz equations (see e.g. §): 

where S is the absolute value of a spin localized on a magnetic ion. In terms of the canonically 
conjugated fields vr and 0, eq. ([16|) can be rewritten in the Hamiltonian form: 

'''''^^^■'^ = s^)' ^"^ 

-RSe,.#>(x,() = -4^KA7r(x,(). (18) 

OTT^X, t) 
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In the harmonic approximation equations ([T7| , [l8D imply the dispersion relation for the 
spin- wave modeEl in a XY magnet: 

6^(k) = A {j^' + 9^^ = (k' +Po^) , (19) 
where c = vX/ is the spin- wave velocity and po = g/J. The out-of-plane anisotropy A 



affects the dynamics in the long- wavelength limit k px = y X/ J considered in this article. 
Thus, px is the upper cut-off momentum in our theory. A spin wave with k ^ pq has the 
phonon-like isotropic spectrum e = ck. The range Pq k <^ px will be called acoustic shell 
(^-shell). At lower momenta p <^ po the spin- wave spectrum is dominated by the dipolar 
interaction: e(k) ^ Cy/pok sin 6, where 6 is the angle between the direction of the spontaneous 
magnetization and the wave-vector. This range of momenta will be called dipolar shell {V- 
shell). The effect of the presence of the dipolar force in the 2D XY-magnet is not limited 
to the change of free spin-wave spectrum. As we have mentioned earlier, it leads to strong 
spin-wave interaction and to a crucial transformation of the spin propagation. 

Without dipolar forces the dynamics of the 2D XY ferromagnet is well described by 
non-interacting spin waves. The high level of fluctuations leads to a strong temperature 
dependence of the dynamic spin correlators, which have algebraic character, just as static 
onesiiii. The renormalized two-point spin correlation function features the pole with the 
temperature dependent power exponent. The dipole force suppresses such strong spin fluc- 
tuations, but not entirely. 

In addition, the dipolar interaction induces decay processes. As a result finite spin-wave 
life-time r(k) or the width of the level 6(k) = r^^(k) appears. In 3D at low temperature 
and at small momentum |k|, the width b(k) is much smaller than ix'(k). In 2D, however, the 
interaction is essential and must be considered seriously. 

To take into account the dissipation induced by thermal fluctuations at a temperature 
T, we introduce a phenomenological dissipation functional^: 

R[(j)]= ydtdWx'i?(x-x')0(x)(/)(x'). (20) 
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Eliminating vr from equations ([T7| , [T8D and adding a proper dissipation term, one obtains a 
following equation for 0(x, t): 

1 



^62(k)0(o;,k) 



where h{'x,t) is the external magnetic field, and the interaction part of the Hamiltonian is 
determined by equation (||). We have introduced the random noise 77(x, t) in equation (^Tj). 
The noise, in effect, generates dissipation. As usual, the random noise is assumed to obey 
the Gaussian statistics. Its correlation function is determined by the fluctuation-dissipation 
theoremii: 

iv^V-^) = 2TR{k) (22) 

Here -R(k) is the Fourier-transform of the function i?(x— x'). The dissipation in the exchange 
ferromagnet vanishes in the long-wavelength limitJil: -R(k) = 6k^. In 2D XY dipole magnet 
the dissipation does not vanish in the long- wavelength limit. We notify: 

fl(k) = (23) 

We emphasize that the finite life-time F is determined self-consistently by the processes 
of the decay and scattering of spin waves. We neglect the spin-wave-electron and spin- 
wave-sound interactions. The first interaction is not weak, but the Fermi-velocity is much 
higher than the spin-wave velocity and the spin-electron interaction is not effective for long 
wavelength. Even the sound velocity can be much larger than the spin-wave velocity, since 
the latter is proportional to small VX. 

Our aim then is to calculate the linear response function G{'x, t) to a weak external field 
/i(x, t) (see equation (|2T|)) averaged over the thermal fluctuations //(x, t). We apply the 
Janssen-De Dominicis functional methodi to reformulate stochastic equation (0) in terms 
of the path integral. The probability distribution for the noise //(x, t) is: 



W[rj\ ~exp 



~j d^dVy dtr/(x, t)R-^{-K - x')r7(x', t) 



(24) 
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Following the standard dynamic field theoryi^'i, upon averaging over the noise distri- 
bution and introducing auxiliary response field 0(x, t), one can reduce the solution of the 
stochastic differential equation (pT]) to the calculation of the dynamical partition function: 



Z[3,j]= /P[0]I?MexpL7[0, 



d^x / dt 



+ 



(25) 



and its derivatives over the currents j and j. Here J'{(j),(l)) is the Janssen-De Dominies 
functional (JDF). 



(26) 



J[^, (p] = J d'x J d'x'J dt0(x, t)TR{^ - x')0(x', t) 
-<^(x', t) 0(x, t)5(x' - x) + (x - x) + i?(x' - x)9,0(x, t) 

By differentiation of the JDF over j and j one can obtain any correlation function. 
We define the inverse bare propagator Gq ^ 



-2TR(k) 



uP _ £f(k) 

A A 



iujR{\i) 



\ 



) 



(27) 



so that the harmonic (" Gaussian" ) part of J can be written in the form: 



1 

'2A 



From eq. (^) one finds the bare propagator: 

/ 

Go = 



(a;,k), </)(^,k) ) Goi(c^,k) 



( ^(_^, _k) ^ 



V 



;-^,-k) 



(28) 



v 



(29) 



^(^^k) 
Go(^,k) Z^o(^,k) 

where we define the bare dynamical response and spin-spin correlation function as follows: 

A 



G'o(^,k) 
^o(^,k) 



cj2-e2(k) -icuA/?(k) 
2TAi?(k) 



(30) 
(31) 



[cu2-e2(k)]'+Cu2A2i?2(k). 

They obey the standard fiuctuation-dissipation relation: Dq = ^2mGo. The same 
relation is correct for the total dynamic correlation D{u!,k.) and the total linear response 
function G{u,k): D = ^ImG. 
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The anharmonic (interaction) part of J is: 



^int=E'^-,k^ = E<^-,k 3 ^ /(k,k2,k3)0k„..0k3. 



.^3 



a;,k 



X 5(k + k2 + k3)5(a; + a;2 + a;3) +4 ^ ^ 5((k, kg, ka, k4)0k2,u;2 0k3,u;3 0k4,, 

a)2,'^3,<^4 k2,k3,k4 



0)4 



X 5(k + ks + ks + k4)5(a; + a;2 + + ^4) , 



We define in a common way the self-energy operator l](a;, k) by the relation: 



(32) 



G-i(a;,p) = Go^(u;,p)-S(a;,p), 



(33) 



where: 



G 



^ G'*(u;,k)^ 



V 



(34) 



(7(0;, k) D(u;,k) 

and G(a;,k) and D(a;,k) are the complete response function and correlator respectively. 
The self-energy S(a;,k) satisfies Dyson equation: 



S(u;, k) = 5 / / G(n, p)A(p, k; 0)G(u; - k - p), 



(35) 



where A(p, k; a;, O) is the full vertex. The bare vertex is given by equation: 
Ao(p,k;a;,Q) = 



^-Do(^^,p) G*(fi,p)^ 



A 



1^ 6-0(0, p) ) 



(36) 



We have denoted /(p, k) = /(p, k, p — k). 

Further we consider a limit i?(k)— >-|-0. According to the FDT, the matrix S(a;, p) must 
have a form: 



S(u;,p) 



^SJmS(a;,p) iS*(a;,p)^ 



(37) 



V p) 

where the self-energy function E(ci;, p) is associated with the complete response function 
G(a;,p) by the same relationship: 



G-^(^,p) = Go-^(.;,p)-1e(c.,p). 



(38) 
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In the Appendix A we show that the full vertex A can be approximated with a good 
accuracy by its bare value /(ti;,p) in the low-frequency range. In this respect our theory 
reminds the Migdal theory of the interacting electron-phonon systemi. In the Migdal theory 
the simplification is due to a narrow scale of the energy shell in which the interaction 
proceeds. In our theory we assume that the frequency of spin fluctuations is small instead. 

We notify the real and the imaginary part of the self-energy term as: S = a'^{uj,p) — 
iujb{uj,p). Thus, the Green function (^) reads: 

G-^{uj,p) = - (p) - (cu, p) + iLub{iu, p) , (39) 

while the spin-spin correlation function is: 

D{u, p) = (40) 

[cj^ — e^(p) — a^(u;, p)] + ^^^^^(ti;, p) 

(we have slightly changed the definitions of G and D and referred the factor V2TX^ to the 

vertex) . 



We employ the reduced temperature t = T/AttJ and the ratio gjy J\ = po/yX/J as 
small parameters. The latter means that the ^-shell is much larger than the P-shell. We 
also use the notation L = \og(\/jX/g). 

The main contribution to the self-energy is given by the one-loop diagrams shown in 
Fig. la, lb. Our theory is valid only if the temperature is small: 



tlog(VJA/^) =tL < 1. (41) 

Under this condition the two- loop corrections are small, and the diagram Fig. lb contributes 
to a negligible change of the spectrum ([T9|) ii. Such neglecting of the two-loop diagrams 
(vertex correction) was a major assumption in the so-called mode-coupling methods0. This 
approximation serves well in the theory of the 3D critical dynamics with the dipole force 
being included. Later we prove this assumption for 2D. Thus, the Dyson equation for our 
problem is as follows: 

E{n, q) = ISA^T / / /2(p, q)D{u, p)G{uj + fi, q - p), (42) 
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The functions b{u, p) and a{u!, p) are even in both argumentscl. The imaginary part of 
the self-energy is odd in u: XmS(r2,q) = —Qb{fl,q). Hence, equation for the dissipation 
function reads: 

b{Q, q) = 9\'T f f f\p, q)D{u, p)D{uj + q - p). (43) 

The integrand in (^31) is positive. Thus, the main contribution to b{Q, q) comes from 
the region, where poles of the two D-functions coincide. The function D{uj,p) has poles 
at uj ^ ^^(p) ill the ^-shell. Following the terminology of the field theory, we call the 
surface u"^ = e^(p) the mass-shell. The self-energy in the ^-shell is small as it is shown in 
the Appendix B and we neglect it. Because the dissipation is small, the D-function can be 
represented as a sum of 5-functions: 

D{u;,p)^J:-^6{Au^), (44) 

where Auj± = uj ± e(p) measures the deviation from the mass-shell. After integrating uj out 
from eq. ( |13| ) with the D-functions from (^), we recover the Fermi Golden Rule for the 
probability of the spin-wave decay and scattering processes. 

Looking for the long wavelength quasi-excitations, we need the self-energy at very small 
momenta q ^ Po, which we denote as Eq. We expect the quasi-excitations to be soft: 
Q ^ cq. Here we restrict the quasi-excitation wave-vector q to be directed almost along the 
magnetization: \qx\ q (arbitrarily directed q are considered in Appendix C). The essential 
contribution to the integral in equation (^21) comes from the internal momentum p being in 
the ^-shell and the internal frequency u = e(p). Performing the integration over u with 
the D-function from equation (0), we find (for derivation see Apendix B and note that we 
neglect the real part of the self-energy in the ^-shell, as it is justified in the Appendix B): 

_ c^plt r c^p^dp VLsm^{2i))di) . 
°" 47r J e\p) n- cq cos i) + ibi ^ ' 

where bi is the dissipation function 6(c<j, p) of a spin-wave inside the ^-shell. Note, that 
TZeTiQ vanishes in the static limit = 0. 
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If cq ^ bi, we make the integral over ifj in equation (^Sj) to find: 



So(x) = cVo^Lcos^xexp(-22x) 



(46) 



where defined by the equation cosx = ^/cq, measures the deviation from the mass shell. 
More generally, we introduce a notation r for the ratio Q/cq (r = cosx if r < 1). Then 



Note that T,o{r) ^ —c^pltjA when r — > oo and |So(r)| < c^plt/i at any r. The self-energy 
So(r) is real for r > 1. 

If q is so small that eg <^ 6i, eq. (|45|) implies the g-independent dissipation constant: 



In this calculation we have used the fact that the dissipation of a spin-wave in the ^-shell 
bi depends only on the angle ip between the direction of magnetization and the spin-wave 
wave-vector p which we prove below. 

Now we need to calculate bi{ip). An unusual feature of our theory is that the dissipation 
process in the ^-shell is mediated by an off-mass-shell virtual spin wave. Indeed, the dis- 
persion relation (|T9D does not allow for decay or merging processes. Alternatively, as we will 
show, the dissipation of a spin-wave in the ^-shell, propagating along the direction specified 
with the angle ip (sin?/; = qx/q), is mediated by an internal virtual spin- wave in (9), with a 
momentum of p <^ po and a frequency of < cp, propagating along the direction very close 
to the y-axis (f'^<^l (where sinif = Px/p), to provide a finite attenuation of this state. The 
integration over uj with one of the D-functions in (|^), taken in the form (^), leads to a 
following equation: 



Since u = e(p + q) — e(q), we conclude that u = cp cos^, where ^ = iIj — ip^ip is the 
angle between the vectors q and p. According to fiuctuation-dissipation theorem: D(ep+q — 




(47) 




(48) 




(49) 
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P) = XmG'(ep+q — CqiP). Invoking the definition of tlie angle x for virtual spin- 

wave, we find that x = ^ i'- Substituting /^(O, q) = (1/9) sin^(2?/;), e^(p) = c^p^ + 
c^PoP sin^ if ^ c^p^ + c^p^pif^ and taking into account that Sq from p6| ) depends only on 
X = '?/', we write: 

. - c^plt f sin^^cosV'dpd*^ 

^^^^^ = ^^"^ y sin^ ^ + p,p^^ + • ^ ^ 

Note that the most dangerous region of integration is the region of small p, such that 
v"^ ~ p/po <^ 1 in (50). In other words, the dissipation of a short wavelength spin-wave. 



propagating in the direction ip, is determined by the scattering on the long wavelength 
virtual spin-wave, with the momentum along y-direction, which lies on a specific distance 
off the mass-shell: uj/cp = cosip. The integration over p in (|50D is confined towards the 
crossover region: p ^ Pc = poVtL. 

To find the anisotropic dissipation of a spin- wave mode in the ^-shell, we plug ^o{ip) 
from equation (^61) into (0). After a change of variables {p,(p) {p,i^), given by formulae 
p = poP^cosi? and (p = psimjj sin.'t} cos"^^"^ {—oo < p < oo and < -i? < vr/2), the 
integration becomes trivial and gives: 

/^.3/4 sin3/^(2^)sin(^/2) 
biW = Pit ' cpo YvZ 'i ' (^1) 

where the direction of the spin- wave is limited to the fundamental quadrant: < ip < n/2, 
and pi = r2(l/4)/4v/2^ ^ 1.31. 

Let us return to the range of very low momenta p ^hi/c. Plugging (|5l|) into (^8[) , one 
finds: 

6o = Pocp,t^''L'>'\ 

where (3q ~ 1.24. The condition cpoM ~ &i defines the crossover wave-vector: Pdm ~ 
PiPot^^^ /L^^'^, between the self-energies (^ and (^8]). The dissipation functions ( ^^ and 
If) represent the self-consistent solution of the Dyson equation (E2|, B3[). 



Finally, we verify that the two- loop correction (see Fig.lc) is negligible. There exists 
several diagrams with different arrangements of G and D functions. On each of the two 
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short loops on the diagram Fig.lc there exists at least one D function but there may be 
two of them. We consider only the most 'dangerous' diagram with each short loop having 
exactly one /^-function (see Fig. 2). Note that the main contribution to the diagram Fig. 2 
comes from regions of internal momenta pi and p2 are restricted to the ^-shell. Inside the 
^-shell the Green and the D-functions have strong singularities on the mass-shell. As it was 
done in the Appendix B we integrate in both short loop the internal frequencies uoi and 002 
and find that only non-static term is non-zero: 

S - el d^P^^^P^ ^/^(Pl,P2)/(Pl,Pl)/(p2,P2) 

J e'^(pi)e^(p2) — eg cos 01 — ibi){^l — eg cos 02 — ^^2) 

^"'(Pi-P^) (52) 
(e(pi) - e(p2) - e(pi - pa))' 

where 0i and 02 is the direction of the momenta pi and p2. We assumed that 0i ~ 02. 
Since the three spin-wave processes is not allowed by the conservation laws, the Green 
function G{uJi — uj2, pi — P2) is off the mass-shell (the corresponding last denominator in 
Eq.(|52l) reads: {pi +P2)(0i — (P2Y +Pocos^(0i) > (if pi,P2 ^ Po- Now we can count the 
momenta powers in Eq. (|52D to verify that the integration is convergent towards the dipolar 
momentum po, and, thus has no logarithm. A simple counting of temperatures shows that 
Eq.(^ ~ t'^/bl ~ t^/^. Hence the two-loop dissipation function is b'^ = bo{l + t^/^/L). 
Similar consideration shows that the function b[{ilj,q) — &i(V'), which represents the two- 
loop corrections for bi, is small in t^^^, and is also small in the ratio po/q. 

Having explicit expressions for the self-energy we can analyze the dispersion relation 
cu^ = e^(p) + S(c(j, p) in the range of small u and p. New results are expected for the region 
p < Pc = PoVtL in which Sq becomes comparable with e^(p). In a range of momentum 



Pdm ^ P ^ Pc and angles ijj <C ypotL/p, we find a new propagating soft mode with the 
dispersion: 



u = cp{p^ + poptp^fypoVtL. (53) 
The dissipation of the soft mode grows to the boundary of the region and becomes of 



the order of its energy at -0 ~ ^p^tL/p or p ~ Pdm- There is no soft mode beyond the 
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indicated range. The spin- wave mode persists at p > potL. In a range p <^ pdm and small 
angles a new diffusion mode occurs with the dispersion: 



-ze'{p)t-'/'L'/yPocpo. 



(54) 



The angular range of the diffusion mode increases with decreasing p and captures the 
entire circle at p < potL. 



At the end of this section we would like to remind that equation (^8|) was obtained for the 
quasi-excitation directed along the y-axis. It can be easily checked that in case of arbitrarily 
directed q one must write sm'^{2ip + 20) instead of sm{2ip) in eq. (^) where sin0 = Qx/g- 
However, in this case the integral becomes singular, so one must treat this equation more 
carefully. We will come back to this question in Appendix C. 



IV. RENORMALIZATION OF THE DIFFUSION MODE. 



In this section we concern ourself with the renormalization of the diffusion mode. As 
we established in the previous section, at wavevectors p < pdm the diffusive dynamics term 
dominates (A — limit) in the harmonic part of JD Functional (27-28). The interaction 
between 'diffusons', given by the anharmonic part of the JDF (32), effectively 'renormalizes' 
the 'diffuson' dispersion at very small wavevectors p < Pa Pdm- We shall determine the 
anomalous diffusion onset wavevector pa in the end of this section. 

To simplify further calculations, we introduce a scale transformation of the fields 
(J^/T2)-V40^ (Jc//T2)-i/^0. In these notations the JDF (28,32) is: 



-k a— 0a;,k - \ ak + — 

\ in \ Ofc 



y\ 



J- 



—w 



\Py\ 



a;,k 



n,p 



Q,p 



where Fq = JF, a = \ J / g and 



w = (T'g/.PY'\ 



(55) 



(56) 



Note, that one can get rid off the spatial anisotropy charge a by rescahng the x coordinate 
only. 

Assuming that all terms in the r.-h.s. of JDF (|55D have the same scaling dimensions, equal 



to zero, and accepting the dimension of /cj^ Aj^ = 1, A° = and A^/p = 0, we immediately 
find from the second and third terms that A°^ = 3/2. Comparing the second and the 
fourth terms, we get A° = 2. The first and the fourth terms give A^ + A° = A^. From the 
condition that the scaling dimension of the first term is zero, we have 2A^ + A°+ 1 + 3/2 = 0. 
Finally, comparing the second and the fifth terms, we obtain A° + A° + A^ + 2 = Now 
it is straightforward to find all bare scaling dimensions: A^^ = 3/2, A° = 2, A° = 1/4, 
A^ = 9/4, A° = 1/4. (Note that the notation A°^ stands for the scaling dimensions of 
0(k, cu) and 0(k, cu), while A°^ denotes the scaling dimensions of 0(x, t) and 0(x, t).) 

According to the standard renormalization group procedure, we introduce renormaliza- 
tion constants: (p = Zl^'^cpR, (j) = Zy'^lpn, a = Zaa^, I/Fq = Zi/tI/Tr and w = Z^wr, 
where 



w 



l^l^w. (57) 



First we note that the fields (f) and (f) have the same renormalization constants, which imme- 
diately follows from the fluctuation-dissipation theorem. Indeed, according to this theorem 

D{x,y) = {(j){x)(f){y)) r^ImG{x,y) =Im{^{x)(t){y)). 

Let us divide the Janssen-De Dominies functional into two parts: 

J = Jr + a J. 

The first one is the " renormalized" functional JTr 



J ~ ,1/4 kxky ' 

-a^— 0a^,k - WrI ' -rrr 

i- R I fty I 



R 



PxP y 

ui,k n,p \Py 
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[01] 






2 









(58) 



The second one I\J contains the counter-terms: 



R 



aniZ^Za - l)kl + {Z^Z'^ - I] 



- aR^{Z^ZaZi/r - l)(j)Ru,,]i - (Zn^zf^ - 1)wrI^^'^^^ 



afi\ky\ 

•3/2 



^Rui,k 



R ^ \k„ 



PxPy 



^Rn,p 

















2 









(59) 



A simple power counting shows that all corrections are only logarithmically divergent 
in the dimension of 5/2 (see section II). Hence, in what follows the e-regularization scheme 
with e = 1/2 is assumed. 

Evaluating the diagram shown in Fig. la up to the second order in ky (using the bare Go 
and Dq functions) one finds the one-loop correction to the second term of the JDF: 



In order to cancel the divergency, we put: 



^^«i?*R/^/^y ^ - {Z^Za - l)an. 



rv ry 18 ~2 ,1/2 I 



h ' 



(60) 



The same procedure for the fourth term gives: 

1 



dk„ 



(61) 



In the Appendix we show that both the three-leg and four-leg vertices do not have one-loop 
corrections. Hence: 



Z,,,Z 



3/2 



(62) 



One can easily see that the one-loop correction to the term in eqn (^) vanishes as 

well fseellHH). It means that 
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Z-^Z^ = 1. (63) 



Equations (^^^^3|) have a following solution: 



Zvr=l + ^="^'"VJ. (65) 
Z.= (l + JL-i'"' / ^) ■ (66) 



2561 " J kl" 



Next we introduce the Gell-Mann-Low /5-function 3 = a and the Callan-Simanzik 

anomalous dimensions 7a = , 7r = f^^R^^g^ ^^id t],/, = z^^§^- We denote by the 

scale at which the coupling constant is equal to wr and denote by A the scale at which the 
coupling constant is equal to w. From ( p^ ) and the definition of Z^, one has: 



./i/ V 1287r k^/\ 
And finally, we find the /5-function, which coincide with that found in the statics cas 

1 27 

/?(*H) = -^*i?, + Y^*| (68) 
The fixed point of tlie renormalization group flow is 



~ * 2 



32n 



27 

After performing this procedure for anomalous dimensions one gets: 





9 


7a = 


—w 

QAtt 




7 ^2 


7i/r= 


327r 






^0 = 


6471 ^ 



(69) 
(70) 
(71) 

Using the result for w}^, it is straightforward to find 7a* = —1/6, 7i/r* = 7/27 and rj^j,* = 1/6. 

The long range limit of the functions G and D is found from the Callan-Symanzik 
equationii: 
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1 ( X t \ 

= (a;2 + |^|8/3)7/8 f [j^^ j (^2) 

1 ~ ( X t \ 

where f{x,y) and f{x,y) are arbitrary functions. 

In the static hmit {t = 0) the exponents in the correlation function (^) are exact, as 
was previously found by one of the authors^: 

D{v) = (0(r)0(O)) ~ y^'/') (74) 



For the Fourier-components of the Green function eq.|T^, we find in the region of anoma- 
lous diffuson: 

G(o;,k) = /i (-^, -^J . (75) 

\ thy Ky / 

The anomalous dispersion of the diffusion mode announced in the Abstract follows from the 
last equation. We see that the static dipole contribution is not renormalized in dynamics, 
as was suggested in (011). We also note that the exchange coupling acquires an anomalous 
dimension Aj = 1/3, whereas the dynamic term uj/T acquires anomalous a dimension 
+2/27. Taking into account the anomalous dimension 71/r = 7/27, we conclude that the 
anomalous dimension of is 7^ = 1/6. The interaction between the diffusons in the scaling 
limit reduces the dissipation or, in other words, hardens the diffuson. 

Now we can estimate the wavevector pa, an upper boundary for anomalous diffusion. We 



assume that temperature is small (^If) . Initially, according to Eq.|56|, the bare vertex wq = 
VTg^^^ / J^^'^, and is also small. Under the renormalization flow, the vertex grows with 
the inverse wavevector as the power 1/4 (pTj). The RG flow starts at pdm and approaches 
the fixed point at the root of the Gell-mann-Low function /3{wr) = 0. Invoking eq.|68|, we 
find the fixed point solution wr ~ 1. Thus, the wavevecor pa is defined as the wavevector 
at which wr ~ 1: 

WoiPDAl/PaY^' ~ 1. (76) 
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We see that 

Pa ~ T'T'/\ (77) 

that is very small. Even if t is not small, Pa ~ PdmWq ~ po{ga/J). It is much smaller than 
Po- 

V. SUSCEPTIBILITIES. 

In this section we find the susceptibility to the magnetic field directed along the average 
magnetization (S) (y-axis), the so-called longitudinal susceptibility. We consider the mag- 
netic field in the form H = Hq + 5if(x, t) where Hq is independent of x and t. When an 
additional magnetic field 6H is imposed, a new vertex 5/itj^k[00]-w,-k emerges in the JDF 
(p6|) (we denote h = gcfiBSH; 6h = qg^ibSSH). It leads to a correction 5D{ijj,'k) to the 
correlation function D(ti;,k) 

5D(xi = X2,ti =t2) = /// fh{uj,q)Do{n,k) [Go(^^ + cu,k + q) + G*(^]-a;,k-q)], 

where DQ{uj,'k) and 6*0(^^5 k) are taken from (pO|,pl]) . 
By definition, the susceptibility x is: 

'^^"^ ^ D(x,=X2,t, = t2). (78) 



2 5/i(a;, k) 
Hence, 

X{u, q) = / / D^in, k)Go(^^ + cu, k + q). (79) 

In the most interesting case, when q = 0, all integrals can be evaluated and the final answer 
is: 

( M , 1^(3/4) ^ ^J^y/^ 2T 



X 
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(80) 



In the limiting case = 0, the susceptibihty reads: 

X = const ■ h~^^^. 

This result has been found earlieiS. 

VI. CONCLUSION. 

In conclusion we discuss how new modes can be observed in experiment. The new modes 
appear on a macroscopic scale of length of the order of magnitude l/xm. Even rather weak 
in-plane anisotropy can suppress or disguise the new modes. Therefore, we can expect that 
the new dynamics will be observed in films with very weak in-plane anisotropy. The best 
known candidates for this role are films grown on hexagonal substrates. The hexagonal 
anisotropy is naturally weaker than the tetragonal one, because they are proportional to 
a higher degree of the relativistic parameter. Besides, the hexagonal anisotropy totally 
vanishes at large distances in a range of temperature from {4:/9)Tbkt till Tbkt, where Tbkt 
is the temperature of Berezinskii-Kosterlitz-Thouless transitionilii. The simplest idea is to 
use the (111) face of FCC crystals, such as Ag, Au, Cu. An iron film on the (111) face of 
Ar has been grown by S.Bader and coworkers^. Recently the Ru film has been grown on 
the hexagonal graphite substratelil. Thus, 2D ferromagnets with exact XY symmetry are 
available. 

The next question is: what dynamic effects can be observed and at what conditions? 
As we have noted earlier, if the in-plane anisotropy is not especially small (less than IK in 
energy scale), the only opportunity is to use a six- fold substrate in the range of temperature, 
not small in comparison with Tbkt- In this situation the observation of propagating soft 
modes described in section III seems to be improbable, since all they require t <^ 1, i.e. 
T ^ Tc. However, the anomalous diffusion can be observed even at T ~ Tc, given a 
sufficiently large scale of length (> lO/xm). The best way to observe it is to apply a short and 
inhomogeneous pulse of magnetic field, and follow when the secondary signal will arrive to 
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fixed indicators. Such a picosecond-pulse technique has been recently used for investigation 
of the film dynamics!!. We propose to use the same pulse technique to different films. 

A quick estimation of the time t^^y needed for the secondary signal to reach the indicator 
at distances and Ly along x and y directions shows: tx,y ~ {h/ Jpoa)(Lx,y/ a)^^'y , where 
a is the lattice constant and A'J. „ are the anomalous-diffusion dimensions for axes x and 
y respectively. Thus, ~ h/ ga^L^/a)'^'^^^^ ~ lO'^'^sec and ty ~ h/ga(Ly/a)'^'^^'^'^ ~ Isec, 
where we have assumed ~ Ly ~ Icm. 

The retardation time for the secondary signal is much longer than the time for the 
primary signal propagation. The strong size and direction dependence of the propagation 
time can be used for detecting of the anomalous diffusion. 
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APPENDIX A: ONE-LOOP VERTICES CORRECTIONS 

First we note that the JDF ( ^5]) has the following nontrivial symmetry^. 

(j){k,uj) ^ (j){k,uj) +e6\k)6{uj), (Al) 
0(k,a;)^0(k,a;), (A2) 
kx ^ kx — eawky, (A3) 

ky^ky. (A4) 

The partition function (^) must have the same symmetry 6Z[j,j] = 0. By a standard 
procedureS, we find the implications of the symmetry ( |A1[ - [A4| ), known as Ward-Takahashi 



identities to the so-called "generating functional for proper vertices" T[(f, 0] (It is the Leg- 
endre transform of In with respect to the fields ip and 0): 
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5T 

aw 



dkx k, —uj 
d 5T 



0. (A5) 



dkx 5<^(— k, —uj) 

Now, writing down Taylor- like expansion of T[(p, (p\ over fields ip and and plugging it 
into (|A5|) , we find: 

{Q,k,-k)-awky^T^^ik,-k) (A6) 

OK 



where k is used for (k, uj). It was proven^ that there is no corrections to the term k1/\ky\ due 
to analyticity. According to ( |A6|) , it means that corrections to the three- leg vertex vanish 
if we put one of the two frequencies, which the vertex depends upon, equal to zero. Hence, 
the corrections must depend on the product of the two tus and are small in the framework 
of Section 111. Considering loop-wise expansion of T in the same spirit, one finds: 

oo 
n=0 

Looking only at the divergent parts of the corrections r'^^-' and T^^\ we see that, because 
r'^'') is just the bare action and has no divergencies at all, the one-loop corrections to the 
w-independent three-leg vertex do not diverge. As a consequence of this fact we obtain 
equation (|62|). 

APPENDIX B: CALCULATION OF SELF-ENERGY 

Let us start with the Dyson Eq.(42) corresponding to a one-loop diagram on Fig. la. 
We may integrate the Eq.(42) over the contour in complex plane uj such that the poles of 
G-function are outside the contour: 

So = ^ / f{v - q/2, P + q/2) [{{^ + e(p - ^/2)f - e^(p + q/2) + ih,)'" 

+ {{Q - e(p - q/2))2 - e^(p + q/2) - ^h,y^ 

^ ^ J e^(p - q/2) 

It is convinient to change in the second G-function p on — p: 
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(Bl) 



pIcH r /2(p-q/2,p + q/2)d2p 



27r J O + e(p - q/2) - e(p + q/2) + ifei 

(B2) 



6-2(p-q/2) ^ f-'(p + q/2) 



{n + e(p + q/2) + e(p - q/2)) - e(p + q/2) - e(p - q/2))_ 

Keeping only the lowest order in small momentum q and frequency Q and after a simple 

expansion in the brackets we find: 

^ ^ pIcH r P{P:P)d^P ^ 

° 27r y O + e(p - q/2) - e(p + q/2) + ibi 

-2n + 4(e(p + q/2) - 6(p-q/2)) ^ ^^^^ 



464(p) 

At this point we separate the above integral into the static Q-independent part and the rest 
'dynamical' self-energy. The static Q-independent self-energy (it always real) reads: 



-'St 



-d'p. (B4) 



27r J e4(p) 

Now let us take into account the static self-energy given by the diagram on Fig. lb: 

^ PocH [ pI-pI ^2 .p... 

Comparing Eq.(B4) (remember that /(p, p) = pIpI/p"^) and Eq.(B5) we conclude that these 
cancel each other in the ^-shell. Thus, we have verified explicitly that the Q, q-independent 
part of the self-energy is strictly zero as guaranteed by the Ward identity due to the rotation 
symmetry of the system (see Appendix A). To get non-zero static self-energy we have to 
expand self-energies like Eqs.(B4,B5) in powers of the transferred momentum q. This will 
be done in the Section IV using statics field-theoretical technic not dynamical as in this 
Appendix. The result is that q-dependent static self-energy only matters when the spin- 
wave momentum q is so small that anomalous diffusion sets up. On shorter wavelengths 
like q < Pdm we may safely neglect the contribution of Eq.(B4,B5). 
Now we return to the Q-dipendent dynamical part of the self energy: 
V/ N P¥'tnc^ f sm\2(/)) dpd<P 

The integral over p gives exactly the logarithmic factor L: 
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E„(x.,)^g^Ln/ „-'"^^7'/ d0 (B7) 

TT J il — q COS + ZOi 

The imaginary part of Eq.(B7) could be easily found in the limit cq ^ hi. In this case 
we use the formula 

P- + m5{x) 



X — iO X 

to find the $J-part of Eq.(46). The K-part is also simple to calculate: 

pIcH fl f sin^ (f) cos (f)[{-fl + q cos (p) +D] 



^^oix, g) = / <Pcosn^-^l + qcos<P)+^ ^^^ 

TT O 7 — COS 6 



q J Q — q cos (f) 

and then simplify: 

pIcH [ sin^ 0COS 



UMX,,) = ^-^LtL n^l^M. (B9) 
TT q J il — q cos 6 



q J Q — q cos 4> 

We repeat the same step: 



TT g"' y il — q cos (f) 

with further simplification: 

3f^So x,g) = L— / -d(p L—— Bll) 

TT q^ J il — qcoscp n q^ 2 

One could easily continue the same procedure to find 

KSo(x,g) =pIcHL^{2^ - 1) =pIcHL cos^xcos{2x), (B12) 

which is exactly the Eq.(46) of the 111 Section. 

Finally let us show that the self-energy in the ^-shell is neghgible. We first assume that 
it does neghgible: a^(a;, p) <S e^(p)- Then we consider the self-energy Eq.(B2), provided the 
external frequency and momentum lies on the mass-shell inside the .A-shell: Q = e(q). Let 



also ^JX/ J ^ q < Pq. One could easily verify that the main contribution comes if g <^ p. In 
this case we may use the Eq.(B6): 

Eo = ^2cg / ""'(''^^ ^ (B13) 

^TJ" J cJq"^ + poq'^/q — qcos(j) + ibi p 
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Integration over p gives \og{JX/ J/q), whereas the integration over the relative direction of 



internal and external spin-waves gives factor yq^/poQx- Thus, the real self-energy 

a^(n, q) ~ tLc^poy/p^. 
At small t this result justifies our neglecting of self-energy in the ^-shell. 

APPENDIX C: THE ANGULAR DEPENDENCE OF THE SELF-ENERGY 

OPERATOR. 



In this appendix we analyze the angular dependence of the self-energy part found in 
Section III. Repeating all arguments, one can find that eq. (|45|) must be written in shghtly 
modified form: 

c^plt f c^p^dp Vt sin^(2?/' + 2</))d?/' 



47r J e*(p) i I — cq cos ip + ibi 
Now again, assuming cq bi, one finds: 



c^pltL 



(2r2 - 1 - 2rVr2 - l) cos(40) + - , ^ sm'^{2^) 
V ^ 2 vr^ — 1 



(C2) 



where r = Q/cq. In the case (j) = and r = cosx the formula (^q ) is recovered. 

The main contribution to the integral in eq. ( ^9|) comes from a region of very small 
angles 0. Hence, the correction to the Sq we have found above is not important in this 
calculations and equation (|50D still holds. 

By the next step we need to plug bi from (^) to the nonzero-angle form of equation 



q, 0) = inc^phL [ 
Jo 

Evaluating this integral, one finds: 



sin2(2V' + 20) 
4:71 ibii^ip) + Q — cq cos ip 



E(|f2| ^ cq) = iVtcp 
cq) = iQcp 



-3/5 



cos^(20) — a 



.Pdm, 



sin2(20) 



-3/5 



cos2(20) -« ( 1 sin2(20) 

\cpdmJ 



(C3) 



(C4) 
(C5) 
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where we denote p = f3ot^/^L^'^pQ, pdm = At^/^L-^/Vo and a = ~ 0.54. In 

what follows we will also use the notation = pQ^/tL. While equation (|C2|) for S holds for 
cq, Q ^ bi{ip), equations (|C^) and ( |C5| ) are valid in the opposite case cq, Q <^ bi{ip). 

Now we can analyze the dispersion relation co'^(q) = e^(q) + S(ci;, q) more accurately. In 
experiment usually t ^ 1 so ^ P ^ Pc ^ Pdm- Easy, but tedious calculations show that 
there can exist up to 9 asymptotic regions in the momentum space with different dispersion 
relations. 



,2 




if Pdm 



siii3/2(2(^) siii(0/2) 



<^ q <^Pq sin 



2 







if PDMiTr/2 - 0)1/2 < g < 
if po0^ -^iq -^ipc 

if PDAf0^/^ < g < ?^0^ 





if - 0)10/3 < g < ^(7r/2 - 0)^/2 

if PDM&'{^/2 - <Pf''' « q « ^0^/^(vr/2 - 0)10/3 



where < < 7r/2 



In figure 3 we show the case t ~ 0.3, L ~ 1 



and Po = 1- 
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Figure Caption 

• FIG.l. a,b) The main contribution to the self-energy. The functions G and D are 
given by (|5^). The three-leg vertices in fig.l.a and the four-leg vertex in fig.l.b are 
from (P^). c) Two- loop correction to the self-energy. Momenta of internal lines are 
indicated. 

• FIG. 2. The most 'dangerous' two-loop diagram. 

• FIG.3. Maps of the regions with different dispersion relations. 
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